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Abstract 

Spin foam models are the path integral counterparts to loop quantized 
canonical theories. In the last few years several spin foam models of 
gravity have been proposed, most of which live on finite simplicial 
lattice spacetime. The lattice truncates the presumably infinite set 
of gravitational degrees of freedom down to a finite set. Models that 
can accomodate an infinite set of degrees of freedom and that are 
independent of any background simplicial structure, or indeed any a 
priori spacetime topology, can be obtained from the lattice models by 
summing them over all lattice spacetimes. Here we show that this 
sum can be realized as the sum over Feynman diagrams of a quantum 
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field theory living on a suitable group manifold, with each Feynman 
diagram defining a particular lattice spacetime. We give an explicit 
formula for the action of the field theory corresponding to any given 
spin foam model in a wide class which includes several gravity mod- 
els. Such a field theory was recently found for a particular gravity 
model Our work generalizes this result as well as Boulatov's and 
Ooguri's models of three and four dimensional topological field theo- 
ries, and ultimately the old matrix models of two dimensional systems 
with dynamical topology. A first version of our result has appeared 
in a companion paper [Q]: here we present a new and more detailed 
derivation based on the connection formulation of the spin foam mod- 
els. 



1 Introduction 

The spin foam formalism p, |, ^, ^ ||, has emerged in the last few 
years as an elegant synthesis of several approaches to quantum gravity and 
diffeomorphism invariant theories more generally. It can be viewed as a 
"path integral" formulation corresponding to the canonical loop quantization 
framework |]T0|, 0, |T2|, |13|, |T^ and also as an extension of the simplicial 
framework for topological field theories |16|, |l^ [l^, which allows more 
general, non-topological, field theories to be represented. 

Spin foams are coloured two dimensional complexes consisting of two di- 
mensional faces (of arbitrary topology), joined on edges of valence > 3. Faces 
are coloured with non-trivial irreducible representations of a "gauge group" 

while edges carry "intertwiners" - G invariant tensors in the product repre- 
sentation formed by the representations on the incident faces. In a spin foam 
model the spin foams are regarded as histories of the physical system and 
assigned quantum amplitudes. In models of gravity a spin foam defines a dis- 



crete spacetime geometry |T5|, |T9| ; Spin foam models of gravity incorporate 



the discreteness of the geometry of space first uncovered [EO, pTl p2|, ESl E| 



in loop quantized canonical theory into a spacetime sum over histories for- 
malism naturally suited to a 4-diffeo invariant theory. This is very appealing 
because the discreteness, which is not "put in by hand" but arises naturally 



^ We will see that G corresponds to the local gauge group or fibre bundle structure 
group commonly refered to as the "gauge group" in discussions of Yang-Mills theory, as 
opposed to the group of all gauge transformations of the fields in spacetime. 
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in this framework,^ promises to remove the ultraviolet divergences found in 
perturbative theories of quantum gravity and matter fields. 

A number of spin foam models of Euclidean^ |2B[ and Lorentzian 



31] quantum gravity in four spacetime dimensions have been pro- 



posed. Three dimensional general relativity [|T^, ^, ^ as well as four dimen- 



sional BF theory [0, |T8|, Q can also be given a spin foam formulation. Even 
hypercubic Yang-Mills theory can be expressed in this framework. In fact 
any lattice model in which a connection with a compact gauge group forms 
the boundary data can be translated into a spin foam model p. 



Of the spin foam models of gravity referred to above all but P and ^0 
are simplicial lattice models. A shortcoming of such models is that their 
predictions depend on the simplicial complex chosen to represent spacetime. 
Topological field theories, such as three dimensional general relativity or BF 
theory, can be formulated on a triangulated manifold in such a way as to 
be independent of the particular triangulation of the manifold chosen, essen- 
tially because such theories have only a finite number of degrees of freedom 
associated to global topological features of the manifold. Four dimensional 
gravity models on the other hand are expected to have infinitely many de- 
grees of freedom in any given spacetime topology, so any finite triangulation 
(in which each simplex carries a finite set of degrees of freedom as in the pro- 
posed models) necessarily truncates and thus misrepresents the gravitational 
field. 

Some sort of continuum limit must be taken. One approach is to sum 
over triangulations of the spacetime manifold. This seems very difficult to 



^ In using the spin foam formalism, or the canonical loop quantization formalism one 
is making an ansatz, so it is possible that the discreteness is not real because the whole 
framework is not the one used by Nature. In lattice spin foam models there is another 
possibility, namely that an infinite renormalization of the bare theory that is necessary 
to go to the continuum limit wipes out the discreteness of the geometry in this limit. 
Finally it should be noted that in a spin foam model of Lorentzian 2+1 general relativity 



developed by Freidel 1 25 geometrical observables do not have an entirely discrete spectrum 
but there is a non-zero minimal length. This is a result of the non-compactness of the 
three dimensional Lorentz group SO{2, 1) and might be a general feature of Lorentzian 
models. 

These models are not Euclidean field theories in the ordinary sense because essentially 
the exponential of i times the action is used to weight histories. We call these models of 
Euclidean quantum gravity because they are proposals for quantizations of classical general 
relativity with metric signature + + H — h. 
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do in three or more dimensions since it is difficult to identify the topology 
of a given simplicial complex. A technically easier, but also more ambitous 
approach would be to sum over all simplicial manifolds^ and thus sum over 
both triangulations and spacetime topology. It is an old dream of gravity 
theorists to make topology dynamical, thus removing yet another element of 
a priori background structure. It would be interesting to sum even models 
formulated on a continous spacetime manifold, such as those found in [H, 



and ||3^, over spacetime topologies. 

In the present paper we will consider another, even wider, summation. 
Simplicial spin foams live on the 2-skeleton of the cellular complex dual to 
the simplicial complex. In other words, the spin foam faces consist of dual 
2-cells, and the spin foam edges of dual 1-cells. The dual 2-skeleton, the 
two dimensional complex formed by these cells (together with the dual 0- 
cells), thus is spacetime as far as spin foams are concerned. We therefore 
propose to sum spin foam models of four spacetime dimensional gravity over 
all "admissible" 2-complexes - 2-complexes that preserve some simple local 
features of dual 2-skeletons of triangulated 4-manifolds. Specifically we will 
require of admissible 2-complexes that their 0-, 1-, and 2- cells are topolog- 
ically points, line segments, and disks respectively. Furthermore we require 
that the combinatorial structure at each 0-cell is that of a dual 2-skeleton to 
a triangulated 4-manifold. Thus each 0-cell has incident on it five endpoints 
of 1-cells and ten corners of 2-cells with each incident 2-cell corner bounded 
by two of the incident 1-cells. This corresponds to the fact that a 4-simplex 
(dual to a 0-cell) has in its boundary five 3-simplices (dual to 1-cells) and 
ten 2-simplices (dual to 2-cells). (See Fig. |l]). In a sufficiently small neigh- 
borhood of a 0-cell the 2-complex is indistinguishable from a dual 2-skeleton 
of a 4-d simplicial complex. Note however that we do not forbid 1-cells and 
2-cells from being incident several times on the same 0-cell, so some fairly 
strange 2-complexes are included in the sum. 

Our main result is that this summation of a spin foam model over 2- 
complex spacetimes may be realized as the sum over Feynman diagrams in 
a perturbation expansion of a quantum field theory: The spin foam model 
defines a (somewhat unusual) quantum field theory on the Cartesian product 

^ Simplicial manifolds are simplicial complexes that are also manifolds. To be a man- 
ifold a complex must satisfy several local requirements the most non-trivial of which is 
that the boundary of the union of the simplicies incident on a point is a sphere. 
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Figure 1: This figure illustrates simplicial complexes and their dual skeletons 
in two and three dimensions, as an aid to the reader contemplating the the 
four dimensional simplicial complexes and their dual 2-skeletons that appear 
in our discussions. On the left a complex of seven triangles (2-simplices) 
along with its dual 1-skeleton is shown. On the right a complex of six 3- 
simplices (sharing a common 1-simplex) and its dual 2-skeleton is shown, 
with the dual 2-cells shaded in. In both cases the dual skeleton is cut off at 
the boundary of the simplicial complex. 
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of four copies of the gauge group G the Feynman diagrams of which are pre- 
cisely the admissible 2-complexes, with the amplitude of each given, modulo 
symmetry factors, by the spin foam sum on that 2-complex. (Comparing 
with a familiar QFT such as scalar theory in Minkowski space we see 
that the spacetime 2-complex plays the role of the graph of a Feynman dia- 
gram, while the spin foams living on the the 2-complex play the roles of the 
values of the momenta on the edges of a Feynman diagram, which are inte- 
grated over to obtain the amplitude of the diagram). We propose to adopt 
the perturbation series of the field theory, with its symmetry factors, as the 
definition of the spin foam model summed over 2-complex spacetimes. 

To make things more concrete let's consider the amplitude for given 
boundary data on the boundary of a spacetime region. This amplitude is 
the sum of the amplitudes of the histories that match the boundary data.0 
In a spin foam model on a fixed 2-complex J representing spacetime the 
boundary is a graph F = dJ and the sum runs over all spin foams living on J 
matching the boundary data - irreducible representations on the edges of F 
and intertwiners on the nodes. Now we wish to sum the amplitude obtained 
by summing over spin foams also over all J bounded by the same, fixed, F,^ 
which is to be considered part of the boundary data. In the field theory 
picture this sum over 2-complexes is the Feynman diagram expansion of the 
expectation value of an observable that encodes the boundary data. The 
field theory formulation provides us with formal functional integral expres- 
sions for this expectation value, and any other quantitiy of interest, which 
are often easier to manipulate than the original sums. Moreover, a regular- 

This amplitude can be thought of as the quantum probabihty amphtude of the bound- 
ary data in the Hartle-Hawking vacuum. Alternatively, if the boundary is divided into 
past and future parts then the amplitude can be interpreted as the transition amplitude 
from past to future data. 

^ In the sum we are contemplating T need not be the boundary of J in quite the 
usual sense. We allow the possibility that parts of V are glued together instead of to 
J. The situation is similar to that of encountered with shrink wrapped foods in which 
part of the plastic covering is stuck to itself rather than the food item. Returning to a 
more theoretical setting, this peculiar feature of the boundaries we admit is physically 
quite natural. If we divide F into a future half r+ and a past half r_ and interpret the 
amplitude of the boundary data as the transition amplitude from the past to the future 
boundary data, then there is no reason to exclude complexes representing "bubble" time 
evolution in which parts of and r_ coincide, and no reason to exclude such complexes 
from a sum over complexes interpolating between r+ and r_. 
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ization scheme that renders well defined the functional integral would also 
provide a definition of the values of the sums. We shall return to this point. 

Our work generalizes the results of Ooguri [|T^] and De Pietri et. al. 
who obtained field theory formulations for two particular models (summed 
over spacetime 2-complexes) , namely SU{2) BF theory [^] and the Barrett- 
Crane model of Euclidean quantum gravity. Many of the underlying 



ideas go back to the work of Boulatov [33|, on three dimensional gravity, and 



to matrix models of two dimensional theories |^ . 

The remainder of the paper is organized as follows. In §H we define 
precisely the class of spin foam models under consideration and then we 
present and verify our main result: an explicit formula for the action of the 
field theory that defines the sum over 2-complexes of any given spin foam 
model in this class. The Turaev- Ooguri- Crane- Yetter model of BF theory 
is discussed as an example. Throughout this section spin foam models are 
formulated as lattice gauge theories on a special type of lattice. In §^ we 
show how these spin foam models, and the field theory generating the sum 
over spacetime 2-complexes, are represented in terms spin foams. The issue 
of divergences, regularization and renormalization is also briefly touched on. 
In the last section we give an argument suggesting that the field theory is 
actually finite for regulated spin foam models (which is not at all obvious 
from the perspective of a sum over spacetime 2-complexes) and outline how 
the theory can be extended to accomodate more general spin foam models. 



2 Field theory formulation of spin foam mod- 
els with dynamical topology 

2.1 Local spin foam models in a lattice connection for- 
mulation 

In the present section we obtain, explicitly, a field theory defining the sum 
over 2-complex spacetimes for a wide class of "local" spin foam models liv- 
ing on the dual 2-skeletons of triangulated 4-manifolds, or more generally 
on admissible 2-complexes. This class includes all Euclidean simplicial four 
dimensional gravity models except that of Iwasaki . 

Instead of working directly with spin foam sums we will use the connec- 
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Figure 2: Atoms and their relation to simplices are illustrated in three di- 
mensions. Panel a) shows an atom and one complete dual 2-cell of which 
the atom has a wedge. Panel b) shows the 3-simplex inside which the atom 
would live were it part of a dual 2-complex of a three dimensional simplicial 
complex. The "four dimensional" atoms that we deal with in the text differs 
from that shown in that it is identical to the portion inside a 4-simplex of the 
2-skeleton dual to a four dimensional simplicial complex. This means that 
the atom has five one dimensional cells and ten wedges, as opposed to the 
four one dimensional cells and six wedges shown here. 



tion formulation of the spin foam models because this provides the easiest 
route to our result. In §^ we will show how the result looks in the spin foam 
formulation and indicate how it may be obtianed within that framework. In 
the connection formulation a history consists of a connection specified by 
elements of G defining parallel transport along "boundary edges" which run 
from the center of a 2-cell of J to the center of one of its bounding l-cells.[] As 
illustrated in Fig. |^ a) these edges cut J into disjoint, topologically identical 

^ Where the centers of these cells are placed does not really matter. In fact the parallel 
transporters could simply be associated with pairs consisting of a 2-cell and a l-ccU in its 
boundary. 
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pieces which we shall call "atoms", because they can be viewed as the fun- 
damental building blocks from which any admissible 2-complex can be built. 
Each atom contains one 0-cell, five one dimensional cells called "spokes" (por- 
tions of 1-cells of J), and ten two dimensional cells called "wedges" (portions 
of 2-cells of J). It is bounded by a graph made of boundary edges that is 
homeomorphic to the 1-skeleton of a 4-simplex, having five 4-valent nodes 
connected in all possible ways by chains consisting of two boundary edges. 
The 4-valent nodes are the centers of incident 1-cells while the 2-valent nodes 
in the chains are the centers of incident 2-cells. (Henceforth when we speak 
of the nodes of the atomic boundaries we shall mean only the 4-valent nodes 
unless the 2-valent nodes are explicitly included). When J is the dual 2- 
skeleton of a four dimensional simplicial complex the boundary edges are 
precisely the edges where the boundaries of the 4-simplices cut the dual 2- 
skeleton, and the atoms that they cut J into are the portions of J inside each 
4-simplex. Atoms are in this sense the dual 2-skeletons of 4-simplices. See 
Fig. |b). 

The models we shall consider have compact G and they are local. Each is 
defined by a "vertex function" V, a gauge invariant function of a connection 
on Fs. This function evaluated on the connection on the boundary of an atom 
gives the quantum mechanical amplitude for that connection. The amplitude 
for the whole connection on all of J is then the product of the amplitudes 
for all of the atoms. 

If we number the 4-valent nodes of F5 from 1 to 5 and let the indices 
i,j, k, ... range over these numbers then we may indicate the oriented edge 
of F5 from node i to node j by lij and the half edge from i to the center of 
lij (which corresponds to a boundary edge in the boundary of an atom) by 
. V is then a function of the parallel transporters gij = hi^[h/]~^ along 
the lij, where hi^ E G is the parallel transporter along CiK In the model the 
amplitude for one atom, x, is thus V{hrci^[hx/]~^) where h^i-' is the parallel 
transporter along the boundary edge mapped to Cj-^, and the amplitude for 
the connection on the whole spacetime is 

w= n v{K^[K/]-'). (1) 

atoms of J 

(Since each atom contains one 0-cell the product can also be viewed as a 
product over 0-cells of J).0 

® A footnote on vector and matrix notation: In (^ h^.,- indicates the whole matrix of 
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If the model is a simplicial approximation to a spacetime manifold field 
theory, as all the gravity models aim to be, then V is an approximation to 
exp(i [Effective action in 4-simplex]) in terms of the connection on the bound- 
ary of the atom in the 4-simplex. The atomic boundary is a graph in the 
boundary of the simplex so the connection on it is a discrete approximation 
to the continuum connection on the boundary of the 4-simplex. 

As is proper for a local theory regions of spacetime (composed of whole 
atoms) communicate only via boundary data - the connection on the bound- 
ary. A model in this connection formalism is nothing other than a lattice 
gauge theory defined on an unusual lattice in order to bring boundary data 
and locality to the fore. For this reason this class of models was called "local 
lattice gauge theory" in [Q. 

2.2 The main result: the field theory that generates 
the sum over topologies 

We can now introduce the field theory that generates the sum over admissible 
2-complexes of the model defined by the vertex function V. It is a real scalar 
field theory on G"^, the Cartesian product of four copies of the gauge group 
manifold, determined by the action 

m = iom-\v[^p], (2) 

where 

^oM = tAt / ^'(^1' h, h,) (3) 
z ■ 4! JG" 

and 

= ^ / d''% V{hi^[h,r^) tP{h')tP{h2'mh')tP{h^)ij{h'). (4) 
5! 

group elements [^;^ Ji<i<j<5 in a sort of abstract index notation. It should be clear from 
the context when such expressions denote a whole matrix or vector or just one element. 
On other occasions we will use boldface letters to denote matrices or vectors, so the matrix 
of hxi-'s would be written as h^: and the second row of this matrix can be written as h.x2 '- 
Again information that is clear from the context will be left out, with h in one context 
refering to all hs in the whole 2-complex while in another only to the four hs associated 
with a given node on an atomic boundary. Finally we shall use the Einstein summation 
convention on gauge group tensor indices, so any repeated index of this type is summed 
over. 
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Jo is the kinetic term. It is quadratic, but it contains no derivatives. AV is a 
5th order, non-local interaction term, with A the coupling constant in which 
we will expand to get the perturbation series. The scalar field ip{hi, /i2, h^,, h^) 
is required to be symmetric in its four arguments. 

A wavefunction 9 of the connection on boundary F is represented by the 
observable]^ 



" ^ ' aSnodes of V 



The arguments hj' of ip are the parallel transporters along the boundary 
edges ej' from the node a toward the four neighboring nodes (indexed by 
6); hr represents the whole connection on F, i.e. the vector of all the par- 
allel transporters along boundary edges in F; and sym{T) is the number of 
symmetries of F, i.e. the number of mappings of the set of nodes of F to 
itself that preserve the adjacency matrix. (Since the identity mapping is a 
symmetry sym{T) > 1). 

The main result of the present paper is the following: 

Theorem 

The formal perturbation series of the expectation value of 0*, 



Vt/j e-^ol'^l+^^l^l e*[^], (6) 



in powers of A is 



\n{J) , 

E I ^^*(hr)^(h), (7) 

JeAr 



sym{J) 



the sum over the set Ap of admissible 2-complexes J bounded by F of the 
overlap of the state 6 with the Hartle-Hawking vacuum of the spin foam 
model on J. n{J) is the number 0-cells in J and sym{J) is the number of 
symmetries of J, i.e. of mappings of J to itself that preserve the combinato- 
rial structure of J.|^ Dip is normalized so that / Dip e"^"''^! = 1. 



^ The observable, mentioned in the introduction, that represents a particular boundary 
connection hp is obtained by taking to be a delta distribution on the gauge equivalence 
class of hr. 

The identity mapping is included in this set so sym[J) is at least 1. 
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Proof: The proof is similar to the derivation the Feynman diagram expansion 
in local field theories in Minkowski space. The order A" term in is 



Vijj e 



(8) 



The functional integral is Gaussian. It's value is zero if the number of factors 
of ip in V^Q* is odd. If the number is even the value of the integral is a sum of 
terms associated with each possible partition of the factors of ip in V"6* into 
pairs. Here we imagine that V^Q* has been written as n factors of V followed 
by G* and that the Vs and G* have been written explicitly as the integrals (Q) 
and ® with the factors of ip in the integrals appearing in some definite order. 
The pairings being summed over are the distinct pairings of the elements in 
the ordered sequence of ips appearing in this explicit expression for V"G*. 
The term corresponding to a given pairing is obtained by replacing each pair 
ip{h.)tlj{h') by the corresponding propagator 



Here h and h' are the sequences of four group element arguments of tp, and 
cr[h'] = ^'o-(4)] is h' reordered according to the permutation a. Note 

that the symmetrization had to be introduced in ([T0|) because the integral 
runs over symmetrized ip only. Thus the integration by parts is justified only 
for the symmetrized derivative.^ 

V"G* is an integral over the connections on the boundaries of n separate 
atoms and on F. When each of the pairs of ips is replaced by a particular 
term in the sum ([T2|) for the propagator the delta distributions reduce the 
integral to one over matching connections on the n atomic boundaries and 

Note also that the delta distributions in (|l2|) are normalized Haar measure deltas. 
The definition (|3|) of Iq is a Haar measure integral and the functional derivatives in ( p^ ) 
are Haar measure functional derivatives. 




(10) 



(11) 



(9) 



Khi,K(i))Kh2,K{2))Kh^^K{z))Khi^K(i))- (12) 



(t£5'4 
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r, where in each term the nodes are matched up according to the pairing of 
the ips and the boundary edges attached to each node are matched according 
to the permutation cr G S'4 of the particular term of ([T2|) being employed. 
(Actually, to make the matching of boundary edges corresponding to a given 
a unambigous we need to adopt a convention to specify the unpermuted order 
of these edges on the atoms and on F. Since boundary edges ej' are indexed 
by pairs of nodes a, b the ordering of the ips, which is equivalent to an ordering 
of the nodes, provides a natural choice for the ordering of the boundary edges 
attached to a particular node a.) 

Glueing together the atoms and F in this way produces an admissible 
2-complex J with boundary F. The corresponding contribution to (|^) is 

^U" j-J dhe*{hr)w{h). (13) 

nl 5! sym[L ) J 

To evaluate (||) we must add up the contributions (|1^) corresponding to 
each glueing, i.e. to each to each pairing of ips with permutations a G 6*4 
associated to each pair. Since (p!3| ) depends only on the structure of J we 
can first sum over glueings that lead to the same J by multiplying (|13D by 
the number of such glueings, and then sum over Js. 

In our expansion the atoms and the nodes in each atomic boundary as 
well as in F are numbered, in an order arising from the way in which the (|^) 
was written explicitly as an integral of a product of ips. When the atoms and 
F are glued together to form a 2-complex J this reference numbering defines 
a numbering of the atoms, atomic boundary nodes and F nodes of J. Con- 
versely, if any admissible, n atom 2-complex bounded by F is equiped with a 
"good" numbering then it defines uniquely a glueing of the atoms and F with 
their reference numbering, and thus a pairing of the sequence of ips used in 
the expansion and a set of associated S'4 permutations. Here a good num- 
bering is an arbitrary numbering of the atoms and of the nodes in each atom 
together with a numbering of the nodes of F such that F has the same adja- 
cency matrix as in the reference numbering. From this we conclude, firstly, 
that our expansion generates all admissible, n atom 2-complexes bounded 
by F, so the sum should run over all such complexes, and secondly, that all 
glueings that lead to the same 2-complex arise from good numberings of the 
atoms and nodes in that 2-complex. 

There are n\{5\)^sym{r) good numberings, but renumbering does not 
necessarily produce a new glueing. A good renumbering of J that produces 
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the same gluing is a symmetry of J, because it defines a mapping of the 
set of cells of J to itself that preserves all incidences, and conversely. There 
are thus n\{b\Y sym{V) / sym{J) glueings that produce the complex J so (|^) 
equals 

{JeAr|n{J)=n} ^y'^^y^ I 

If the number of factors of if) in (|^) is odd then n atoms and F together have 
an odd number of nodes and so cannot be matched up to form a 2-complex 
bounded by F. (|Tip therefore covers this CcLSG clS well, because the integral 
(H) is zero while (0) is zero because {J G Ar|n(J) = n} is empty. The 
theorem is thus established. I 

Note that we could have restricted ^4) to be invariant under the 

gauge transformation 

K ^ gK (15) 

without affecting the result. If one writes ip = ip+Atp where ip = jQdg ilj{gh.), 
the gauge average of ip, then one finds /q ['?/'] = IqH'] + -^of^"^] while V[ip] = 
V[V'] and Q[ip] = Q[ijj] because V and 9 are gauge invariant. It follows that 
the integral over Atp cancels between numerator and denominator in the 
expectation value 

(6*) = [Vijj e-^°[^]+^^['^le*[V']/ / Vijj e-^°['^] (16) 



and therefore this expectation value is unchanged if is replaced by ip, i.e. 
is restricted to functions invariant under (^). (In (^ the denominator of 
(|T6|) was set to 1. We may simultaneously set fVxp e~^°^^^ = 1 and leave 
out the denominator also in the ip formulation). In the perturbation series 
of the ip field theory the propagator (|T2p is replaced by one that requires the 
connections to match only up to gauge. However, the gauge invariance of 
V and 6 imply that the integrals over the connection can be evaluated in a 
gauge in which the connections on the atomic boundaries and F really do 
match. The integrals over the new gauge degrees of freedom at nodes then 
simply contribute factors of 1 because the normalized Haar measure is being 
used. 
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2.3 An example: The Turaev-Ooguri-Crane-Yetter model 
of BF theory 

A very simple gauge theory to which we can apply our scheme is BF theory. 
In this theory the action for the boundary value problem in which the con- 
nection is fixed on the boundary is J tr[B A F] where F is the curvature of 
the connection, 5 is a 2-form which takes values in the Lie algebra of the 
gauge group G, and the trace is taken in this Lie algebra. Extremization of 
this action with respect to B requires F = 0, that is, a fiat connection and 
in particular a fiat connection on the boundary. Indeed a naive path integral 
quantization of this problem gives as the amplitude for the boundary con- 
nection a delta distribution with support on fiat connections. It is therefore 
natural to discretize BF theory as a simplicial local lattice gauge theory in 
which the amplitude for the connection on the boundary of each simplex is 
a delta distribution with support on fiat connections or, more concretely, V 
is a delta distribution with support on fiat connections on T^. Here fiatness 
on Fs means that the holonomy around any closed loop of F5 is trivial. The 
only such connections are the trivial connection and its gauge transforms, for 
which 

h^' = m'' (17) 

with Pi ^ G characterizing the gauge at the 4-valent node i while q^^ = g-^* G 
G does so at the bivalent node in the chain connecting i and j. Thus 

VBFihi[h/]-') = Yl|dPkl[5ihAh/r\p^pJ') (18) 

k i<j 

= n/^Pfc n / ^9'™ U^h^^^p^q^^). (19) 

k l<m i^j 

Substituting Vbf into the formula (^ for the interaction term in the field 
action we find 

■■■^ih') (20) 

i^ip^q'') (21) 
(22) 



VBFm = ^l^^^d^°hVBFih,^[h/]-')^Pih')-. 

= ^n/^P^n Jdq'^^{p,q'n-...- 

k l<m 

= jdq'^^l,{q'^)-...-,l,{q'>^). 

' l<m 
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/o[V^] — \Vbf[4'] is precisely Ooguri's [0 action for generating BF theory 
summed over 2-complexes. 



3 Spin foam formulation 

3.1 From local lattice gauge theories to spin foams 

The results of the previous section can also be understood within the spin 
foam framework. This has the advantage, among others, of replacing the 
functional integrals of by more concrete multiple integrals over infinite 
sequences of real numbers, suggesting avenues toward regulating these inte- 
grals. 

Any local lattice gauge theory with compact gauge group G can be given 
a spin foam formulation. The spin foam formulation is conjugate to the 
connection formulation in the sense that if a connection is a history of the 
configuration variables of the model a spin foam is essentially a history of 
the momentum variables - the states in the spin foam picture are a sort of 
Fourier transform of the states in the connection picture. 

The spin foam sum for the amplitude 

J dhe*{hr)w{h), (23) 

of a boundary state ^ on a fixed admissible 2-complex is obtained by expand- 
ing V^(h^) for each atom and ^(hr) on a basis of "spin network functions" 
^, and then integrating out the connection in each term of the resulting 
expansion of 6*w. Each non-zero term in this sum for (p3D is associated to a 
spin foam. The translation of integrals over spacetime connections to sums 
over spin foams is described in detail in 0. Here we only outline its main 
features. They key necessary element for any reformulation of this type is an 
orthonormal basis of the space of distributions on the gauge group manifold. 
Since G is assumed to be compact the Peter- Weyl theorem assures us that 
the matrix elements of the unitary irreducible representations,^ U^^''m,'^{g) 
(with g ranging over G, r labeling the representation, and (m, n) the matrix 

Only one irreducible representation is selected from each class of unitarily equivalent 
representations. The basis of functions on G consists of the matrix elements of these 
selected representations. 
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elements) form an orthogonal basis of distributions. Indeed, since 

with dim r the dimensionality of the representation r, the set {\/ dim rU'^^^ m^} 
is an orthonormal basis. 

Any function of the parallel propagators along the edges of a graph may 
thus be expanded in terms of tensor products of representation matrices of 
the parallel propagators. Notice that a representation matrix of the parallel 
transporter along an edge is a two point tensor, with one index living at each 
end of the edge. The subspace of functions that are gauge invariant at a 
given vertex is thus spanned by functions obtained from the tensor product 
of representation matrices by contracting the indices at that vertex with an 
invariant tensor, also called an "intertwiner" . (An example of an intertwiner 
when G = SU{2) is the antisymmetric tensor eijk, which corresponds to a 
trivalent vertex with each incident edge carrying the spin 1 representation). 
An orthonormal basis of fully gauge invariant functions can be found by first 
selecting at each vertex an orthonormal basis {VT/}/ of the space of intertwin- 
ers for each set r of incident representations that admits intertwiners^ and 
then forming the basis functions by contracting factors ^fdmirl U^^''^rn"{ge) 
for the edges e with intertwiners belonging to the chosen intertwiner bases 
at the nodes. The data labeling a basis function, namely the representation 
on each edge and the basis intertwiner at each node, define a "spin network" 
on the graph,|^ and the basis functions are called "spin network functions" . 

Now consider expanding V^(g) into a sum of spin network basis functions. 

V{g) = Y: A(-)i [n v/^m^]0["-^(g). (25) 

r..,I. i<j 

(Recall Qij = hi^[hj^]~^). 4'i^"\g) is the (normalized) basis spin network func- 
tion on Fs with intertwiners Wj^ at the 4-valent nodes and representations 

That is to say, each r = [ri, r2, r„] such that ri (g) ... (g) contams the trivial 
representation, so that non-zero invariant tensors exist. 

Strictly the spin network consists of the subgraph carrying carrying non-trivial repre- 
sentations - edges that carry trivial representations are left off - together with the "colour- 
ing" data consisting of the non-trivial representations and the intertwiners. In our present 
context of lattice gauge theory this point is of no importance but in the continuum in 
which there is no convenient graph of all possible edges it is essential. 
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rij on the edges connecting the 4-valent nodes. 



i<j 



r.. denotes the whole matrix of representations r^j and I. similarly denotes 
the sequence of basis intertwiners !{. The factor [Yli^j \J dim r^] in (^) has 

been separated from the coefficient A^^--^^- for later convenience. 

Each spin network on the boundary of an atom defines a spin foam on the 
body of the atom: The intertwiner /j on the 4-valent node i is the spin foam 
intertwiner on the spoke connecting the node with the 0-cell at the center 
of the atom; The representation on the chain connecting nodes i and j 
defines the representation on the wedge of the atom bounded by the chain. 



The orthogonality relation (^J) and the orthogonality of the basis inter- 
twiners then ensure that when 6*w is integrated over the connection the only 
terms that remain are those in which these atomic spin foams match up, in 
the sense that the intertwiners on the two halves of a 1-cell of J agree and the 
representations on all the wedges of a 2-cell agree. Moreover, they ensure 
that the intertwiners and representations match at the boundary F where 
the spin foam meets the spin network in the expansion of 6. 

The weight, or amplitude, of each spin foam S in the sum for ( ^3] ) that 
remains once the connection is integrated out isQ 

4^ n n dtmr-y (27) 

O-cells of J j/S 2-cells of J 



They must agree in the sense that if two wedges of a 2-cell have matching ( "coherent" ) 
orientations then they must carry the same representation, whereas if they have opposite 
orientations they must have complex conjugate representations. These representations will 
be discussed more a little further on. 

This formula should be viewed as somewhat schematic some details have been ignored. 
The matching of basis intertwiners on coincident nodes is not between basis intertwiners 
in the same invariant tensor space but rather in dual spaces. If we use a raised index I 
to label the dual basis intertwiners we see that the matching of intertwiner indices A in 
the product ( p7| ) should always be between a downstairs and an upstairs index. For many 
groups one may choose the intertwiner bases so that raising or lowering the intertwiner 
index is trivial, leaving A unchanged, but we are not sure that this is always true. Sign 
factors, not included in our formula, are also present for "psuedo-real" representations that 
some groups have which depend on the relative orientations of the coincident boundary 
edges on which neighboring atoms meet. The proper resolution of these issues is explained 
in outline at the end of this section. 
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where css proportional to the coefficient of the spin network dS (the bound- 
ary of the spin foam) in the expansion of 6, being obtained from the coefficient 
by dividing by a factor of y/dimr for each edge of T. 



3.2 Divergences, regular izat ion, and renormalization 

Spin foam models can have divergences when J has closed 2-surfaces. For 
instance, in the model of BF theory that we have discussed the weight w of a 
connection has redundant delta distributions, and is thus infinite, whenever 
J contains topological 2-spheres. Other closed surfaces can also contribute 
singular factors^] 

In the spin foam sum the infinities take the form of a divergence of the 
sum over representations and can thus be regulated by somehow cutting off 
this sum. This can be accomplished by replacing G with the quantum group 
Gq with g 7^ 1 a root of unity, which has a finite set of inequivalent irreducible 
representations. The spin foams are coloured with the representations and 
intertwiners of Gg, and A is replaced by its natural generalization to such 
representations ||16[ yielding a finite sum. One can also simply cut off the 



sum, leaving the model otherwise unchanged. 

In lattice models of topological theories, or Yang-Mills theories for that 
matter, the continuum limit theory is usually defined not by a sum over lat- 
tices but by the limiting values of the observables of the theory as the lattice 
is refined. In the case of simplicial models the existence of this limit re- 
quires that the model be independent of the triangulation used (at least 
so long as the triangulation is very fine, i.e. has very many simplices). 
This requires that the regulated model be renormalized, because the number 
of divergences, and thus the number of the large, regulator dependent fac- 
tors in the transition amplitudes, depends on the simplicial complex used. 
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The singular factor associated with a closed surface in BF theory is 



"^^{dimr)^ + (— dimr)-^ + k ^^(dimr)^ (28) 



reR r£PR reC 



where x is the Euler characteristic of the surface, k = 1 if the surface is orientable and if 
it is not, R is the set of (unitary equivalence classes of) real irreducible representations, PR 



is the set of psuedo-real irreps, and C is the set of complex irreps (for definitions see |35 ). 
When G = U{1) we obtain a singular factor for every orientable closed surface, whereas 
for G = SU (2) only the four surfaces with Euler characteristic x ^ create divergences. 
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For BF theory on a triangulated manifold this means dividing out a factor 
^(O)reguiated = [Sr ("^^"^ '^)^]reguiated for each independent 2-sphere in J.[^ This 
is the origin of the factor ^(0)^g~i^g(j (with Nq the number of 0-simplices and 
A^^i the number of l-simplices in the triangulation) in the weight of a history 
in simplicial BF theory [|18|. Renormalization is more subtle in the context of 



a general admissible 2-complex because the spheres do not necessarily span 
all the closed 2-surfaces (for instance one can make a 2-complex having a 
torus as it's only closed surface). As far as we know the renormalization of 
the model has not been carried out in this wider context. 

In the present paper we are of course advocating a different approach to 
the continuum limit, in which one sums over 2-complex spacetimes. We have 
not studied the issue of renormalization in this approach. 



3.3 The field theory in terms of spin foams 

Now let's see how the field theory (D) can be expressed in terms of spin foams 
and networks. While this is qualitatively quite straightforward it turns out 
to be a little tricky in detail. The difficulty has to do with the fact that 
the spin network basis is not quite unique. Some conventions have to be 
introduced to define a particular basis. When translating from an integral 
over connections to a sum over spin foams on a particular, given, 2-complex 
one can choose the spin network bases on the atomic boundaries so as to 
simplify the calculations on that particular 2-complex. On the other hand, 
in order to express our field theory in spin foam terms, which boils down to 
expressing it in terms of the coefficients in spin network expansions of V, 6, 
and ip, we have to choose a particular spin network basis on the prototypical 
atomic boundary Fs. This is then the basis on every atomic boundary in 
the complexes generated by the field theory. The resulting infiexibility in 
the choice of basis complicates the calculation of the amplitudes of Feynman 
diagrams. 

We have taken a compromise route. We present the field theory in spin 
foam language using a basis that makes the presentation as simple and sy- 

The expansion on the basis of matrix elements of the Haar measure deha distribution 
on G is 

(5(g) = ^ dim r trC/('^). (29) 
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metrical as possible, obtaining a very neat result. Then we outline how the 
amplitudes of the Feynman diagrams can be calculated by exploiting the fact 
that one is free to change to adapted spin network bases independently for 
each Feynman diagram, making the calculation as simple as in a spin foam 
model on a fixed 2-complex spacetime. 

It may well be that the results could be presented more cleanly using 
a slightly more abstract "basis independent" approach to the spin network 
basis, in which the conventional choices that fix a particular basis are not 
made. Here we will stick with definite bases. 

Let us begin the translation of the field theory by expanding using the 
Peter- Weyl theorem: 

V^(h) =^6W-„nv'fcri f/(^')™,"'(/it). (30) 

r t=l 

The symmetry of if) under permutations of its arguments requires the coeffi- 
cients h to be similarly symmetric: 

6W-„ = 6Hr])-H^^^j VaG^4. (31) 

The free action is easily expressed in terms of these coefficients if we write 
it as /o[V'] = 2^ /g* '^*'^ This is valid because ip is real. We find, using 
(^, that 

^oM = ^E[&^^^"^n]*M^)%. (32) 
The corresponding propagator is|^ 

^^(rOm,^j^(r.)n..^^]*^^ = E ^r,.[r.] (33) 

With 5l = Y{U5l\. 

To express the interaction term V in terms of the coefficients h in as clean 
a way as possible we shall introduce a further representation theoretic tool. 
If ^712 £ G defines parallel transport on the edge from node 1 to node 2 then 
921 = defines parallel transport along the inverse edge from 2 to 1. A 

It can be evaluated by translating the field propagator (■(/'(h)^(h'))o or directly from 
(|3^). In carrying out the Gaussian integral in the second approach the reality conditions 
satisfied by b (spelled out further on in the text) must be taken into account. 
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unitary representation matrix t^mi™^(5'i2) can be expressed in terms of g2i 
as [f/^'"H5'2i)]mi"^ = Wml""^ i92i)]* ■ It is just the parallel propagator from 2 
to 1 in the complex conjugate of the representation r. Thus a spin network 
function of the parallel transporters along the oriented edges of a graph can 
be expressed just as well in terms of parallel transporters of the edges with 
different orientations, provided the representations on the reversed edges are 
replaced by their complex conjugates. 

Now recall that in a spin network basis the representations f/^*"^ that 
may be placed on edges are particular representatives chosen one out of each 
unitary equivalence class. Moreover, in order to keep the definition of the spin 
network basis simple we shall use the same set of representatives on all edges. 
In general the conjugate representation [L/*-^-*]* will not be the representative 
U^''^ of its own equivalence class, it will only be unitarily equivalent to it. 
For instance, all representations of SU{2) are equivalent to their conjugates, 
but unitary irreducible representations of non-integer spin are not real and 
so are not equal to their conjugates. 

We therefore introduce e^^, a unitary matrix such that 

and let = [e^^]*, which is the inverse of in the sense that e^^.e*^''^"^ = 

5^. A few notes on e'^'^'^: 

• It follows from ( p4|) that is an invariant tensor (an intertwiner) with 
m a covector index of the representation r* and n a covector index of the 
representation r. 

• (^) determines e*-**^ up to phase. Consequently £:('■*) equals up 
to a phase factor, and thus when r = r*, e^^^ is either purely symmetric or 
antisymmetric.|^ When r ^ r* one is of course free to choose the representa- 
tive of the equivalence class of the conjugate representation as one likes. One 
could simply choose the conjugate representation itself and use = 5mn- 
We are therefore free to set e*^^-* = e'-'"*-' for all r.0 • e'^'^'^ can be used as a 

For example consider the spin 1/2 and spin 1 representations of SU{2) in the standard 

(-) 

eigenbasis of the generator J^. In the spin 1/2 representation Smn — Cmn, while in the 
spin 1 representation Smli = 5m,-n, so it is antisymmetric for spin 1/2 and symmetric for 
spin 1. Indeed one finds that Smn is symmetric for all integer j and antisymmetric for all 
half odd integer j. 

When r = r* and e^'') is symmetric the representation r is real; If r = r* and e'^''^ is 
antisymmetric r is psuedo-real; Finally, if r 7^ r* r is complex. See p5|. 
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(possibly antisymmetric) metric to turn a representation r vector (carrying 
an upstairs index) into a representation r* covector (carrying a downstairs 
index) according a^*-* = e^^^a^^^^. The inverse, index raising, operation is 
Q^ir)x _ ^{r)xm^{r*) ^ Q£ QQ^j-gg complcx coujugatiou also turns unitarily trans- 
forming r vectors into r* covectors and vice versa. The index positioning in 
all our equations is consistent with this fact. 

The relation (|3^) can be used to define a modified, more symmetric spin 
network expansion of V which depends only minimally on the orientations 
chosen for the edges of V^. Notice that 

uti^'ih.h,') = us^{h)[u!:^%h,)r (35) 

= uS''{h)e^:^U^'Jy{h)e^^^-"^\ (36) 
Thus adopting the notation rij = r*- for i > j we may write 



+3 «<i 



(37) 

The lower indices of the intertwiners appearing in (^) have been raised us- 
ing e. The invariance of e implies that the resulting tensors W^™* are also 
invariant and thus intertwiners for the incident representations r^j. It is easy 
to verify that they form an orthonormal basis of such intertwiners. (|37D is 
thus a spin network expansion of V with the bivalent nodes assigned the 
intertwiners e^xljlji (which is ^^ dim times the unique (mod phase) nor- 
malized bivalent intertwiner). In ( |37D the orientations of the edges manifest 
themselves only in the ordering of the indices in these es. 

Evaluating V is now straightforward. If we use the reality of ijj to replace 
ip by ip* in (Q) the integrals can be carried out using (0) and we obtain 

V = ^ E n xj* w^- n^Ar (38) 

r..,I. k i<j 

V is even simpler when expressed in terms of the spin network coefficients 
in an expansion of the "gauge invariant" field ip{h) = jQip{gh) dg. A spin 
network type expansion of this field can be obtained by inserting into the 
expansion (|30|) of ip the projector 



/ n f/^'*L,"*(^7)f^^7 = E[^i'^"'™'"'"1*W^?"'^"'"'"'"' (39) 
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onto invariant ri r2 ® ra ® T4 tensors. Thus 



^(h) = Ec^'^'nW^i'^'" n ^dimn U^^'^mriht) (40) 

r,/ t=l 

with c'-^'-'^n = &^'"^"^n[W^/''''™]*- III terms of the coefficients c the interaction 
term is 

v=^E^^^-^^-n[c^"'^V]* n^^.. (41) 

r,.,I. k i<j 

It's just A times the result of contracting together the indices of five c* in 
the pattern of a 4-simplex. More precisely it can be obtained from ( P?] ) by 
replacing the intertwiners Wj by [c^]* and setting all the parallel transporters 
to 1. This prescription can also be applied to 6 to express G in terms of the 
coefficiants c*. 

e also allows us to express the reality conditions satisfied by the coeffi- 
cients b: ip* = ip and (|33) implies that 



= [6{r)x^]*e(r)mx^W (42) 

(with e^'")'"" = nti 

We are now in a position to rewrite the functional integral (^) as an 
ordinary integral over the coefficients b (which are finite in number if the spin 
foam model has been regulated by cutting off the sum over representations). 
Moreover we may evaluate the Feynman diagrams of the perturbation series 
in A using the propagator ( p3D and the reality conditions (^21) . 



What we have given is just about the simplest, most symmetric state- 
ment of the spin foam formulation of the theory. However, it is not the most 
convenient form for actually evaluating the Feynman diagrams. In fact it is 
best to first select the Feynman diagram one wishes to evaluate, correspond- 
ing to a particular 2-complex J, and then choose the spin network basis for 
expanding 0* and the factor of V corresponding to each atom. To define this 
adapted set of spin network bases we choose an orientation for each 1-cell and 
2-cell of J. The orientations on the 2-cells induce an orientation on each edge 
of F and each atomic boundary edge that matches the sense of circulation of 
the 2-cell that it cuts. The orientations on the 1-cells define a sign on each 
node, which is positive if the 1-cell is outgoing at the node and negative if it 
is incoming. Now we choose spin network basis functions corresponding to 
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these orientations of the edges. We obtain an expansion of V, rather hke that 
given in (^SD with coherent orientations on the chain of two edges connect- 
ing a pair of nodes. However, we do not necessarily have two outgoing and 
two incoming edges at each node as in (pSj). The intertwiner bases are also 
adapted. At each positive node a basis {VT/^^^} with indices suitably lowered 
for incoming edges, is used. At the corresponding negative node (the neg- 
ative node that lies on the same 1-cell) the intertwiners [lyj'^^'"^'*]* are used. 
Here a G 5*4 is the mapping of the incident edges at the negative node to the 
corresponding incident edges at the positive node. These complex conjugate 
and permuted intertwiners are still orthonormal and have the index positions 
compatible with our convention for the orientations of the edges. 

With these conventions (and with the correspondingly adjusted coeffi- 
cients A and cqs) the amplitude of the Feynman diagram is relatively easy to 
obtain. The result is as described in (pTj): It is a sum over histories consisting 
of representations and basis intertwiners, with only histories in which every 
2-cell y carries a single representation Vy - the common representation on the 
atomic boundary edges that cut the 2-cell and on the edges of T that bound 
it - and each 1-cell carries a single intertwiner basis index /, the common 
value of the indices at the positive and negative nodes on the 1-cell. Such 
a history is a spin foam. The amplitude for each history is the product of 
the coefficients A for the atoms times Cgg and a factor dim Vy for each 2-cell. 
This last factor results from the factors of Exi-xji that appear in (^) and (^), 
associated to each wedge. When the orientations of the edges are changed 
to match those of the 2-cells these are replaced by Kronecker deltas. Then, 
when the propagators for b are substituted in, the Kronecker deltas associ- 
ated to wedges in a given 2-cell end up being contracted in a chain around 
around the 2-cell, thus contributing a factor dimry. 



4 Some closing remarks 

4.1 On the possible finiteness of the field theory cor- 
responding to regulated spin foam models 

The number of admissible 2-complexes increases very rapidly with the num- 
ber of atoms. In fact the number of complexes of n atoms with a given 
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boundary T having m 4-valent nodes is given approximately by 

_ A^-Y5n+m)/2 (5n + m)! 1 1 

iJ,,otHi'HJ)=n} ^ ([5n + m]/2)!n!(5!)"sym(r)- ^ ^ 

This increases as (n!)^/^ for large n, so unless the vertex function V is very 
special indeed the radius of convergence of the series (|^) is zero. 

What prospect is there of making sense of the functional integral (P)? In 
fact it is not unreasonable to hope that once the spin foam model is regulated 
by cutting off the sums over representations (as we discussed in the context 
of BF theory) that the integral (P) can be assigned a finite value by analytic 
continuation. We will suppose that the finite set of representations summed 
over in the cut off model includes the representation r* whenever it includes 
r. Then ij) in the cut off model is a real function on determined by a finite 
set of parameters. These paramenters can be the fe^^'^'^n, but we will use a 
set of linear combinations {xp}^]^ of these chosen so that they are real and 
IqI^j] = x^x. Then (H) becomes 

j ^ d^x e"^"^^ gAVfi- •f5xp,.....xp, e*?!...?™^,^^ . . 3,^^^ (44) 
a function of A which we shall call F{X). For simplicity we will suppose that 



V{^) is real, as is the case in BF theory, and in the models of and [p7 
then also V"^'""^ is real and we see that the integral is convergent when A 
is pure imaginary - the integral is a Gaussian times a function of modulus 
0(1x1™). Can the result be analytically continued away from the imaginary 
axis? 

To see that the answer is quite possibly "yes" consider the simplest analog 
of the integral (H): 

/(A) = r dxe-"'e^"'. (45) 



The coefficients in the formal power series expansion of / about A = also 
diverge as (ra!)'^/^ for large n and, like (|^) the integral converges for purely 
imaginary A. Moreover, / can be continued to the entire complex plane, with 
a fivefold branch cut extending from to cxo. To show this one first writes / 
as /(A) = h[X) + h{—\) where 

h{\) = / dxe-"" e^^ . (46) 
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This integral is convergent when 3?A < 0. It may be reexpressed in terms of 
\~^/^ and the rescaled integration variable y = X^^^x. From this form one 
obtains a series for h in powers of A-^/^ that is convergent for all finite values 
of this parameter, and the claimed result follows. 

The first step in this argument can be repeated for the regulated func- 
tional integral (^41). Let M be the union of rays in on which V"'^'""''^Xni ■ 
■■■ ■ a^ns < 0. Then 

F(A) = H{X) + H{-X) (47) 

with 

H{X) = [ d^x e-"^"e^^(") e*(x). (48) 
Jm 

Again we find that the integral converges for 3ftA < 0. The question is now 
whether it can be continued beyond this domain. It seems likely that this 
would be the case for generic models, especially as it is finite on 9ftA = 0. If it 
is continuable then -F(A) will be defined and finite on an open region in the 
complex plane, possibly including A = 1 which is the value corresponding 
most closely to a simple sum of the spin foam model over all admissible 
2- complexes. 



4.2 Generalizations of the formalism 

Our formalism can easily be generalized to a wider class of 2-complexes. 
We have allowed only atoms which are dual 2-skeletons of 4-simplices, which 
means that they have five spokes each of which is four valent. The field theory 
can be extended so that it generates 2-complexes including atoms with any 
given number p of 4- valent spokes (which are dual to 4-polyhedra bounded by 
p 3-simplices) . All that is necessary is to include a suitable interaction term 
in the action (Q) formed from a vertex function for the new type of atom like 
V was formed from the vertex function V and from the boundary state 6. 
The expression ( P3[ ) for the amplitude of the state 9 on the boundary F of J 
can be viewed as the partition function for a closed 2-complex consisting of 
J and one copy of a new type of atom with boundary — F (F with reversed 
orientation) and vertex function 6*. Adding B* to the action would generate 
all complexes consisting of both this new type of atom and the original type. 

Even with the above generalization all atoms have spokes of valence four 
and thus all bounding graphs of atoms and of complexes have non-trivial 
nodes of valence four only. The theory may be generalized to acomodate other 
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valences q by introducing additional fields depending on q group elements 
instead of four like tp does, and adding a quadratic free action to the total 
action for each such field. 

Any model having a finite set of types of atoms can be handled by these 
means, including of course models with spacetimes of any dimensionality - 
to the extent that this dimensionality can be captured in the combinatorial 
structure of the atoms. 



4.3 Matrix models 

When our formalism is applied to two dimensional models with atoms dual to 
triangles we find that the analog of c^"")^!! is c^^'^^nma = ^(''i-^a) ^dlmvi '^'-i^2*^^iU 

(ri has to equal because ri ® T2 contains the trivial representation only in 
this case, imr is the only normalized bivalent intertwiner) . c'-'^-' is es- 



sentially the matrix of matrix models ||3J] : The triangular matrix model is re- 



covered by choosing the vertex function to be V{gij) = A{diTnr)^ /3\tr U^^\gi2g23g3i] 
so that only one representation r appears in its spin network expansion. Then 
V = A/3!c(:)"^cM"3cM-i. 
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